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Abstract.- The two-dimensional Potts Model with seven states under external field 



(N 

> , 

— , is studied using a cluster algorithm. Cluster size distribution and the fluctuations in 

' the average cluster size provide helpful information on the order of phase transitions. 

(N 

lO ■ Potts model is known to have a very rich critical behaviour and considered as a 

^ . testing ground for both analytical and numerical methods. Recently, the problem 

, is mainly focused on determining the nature of the phase transition, especially the 

weak first-order transition occuring in this model 0, 0] . In two dimensions, the Potts 
model displays a second-order phase transition for the number of states g = 2, 3 
• and 4, and first-order transition for g > 5, where transitions become stronger as q 

increases @, 0. One of the most reliable methods to study first-order transitions is to 
^ ! observe the doubly-peaked probability distribution for energy @, |3] • At the critical 

\ point, the order parameter and the cluster size distribution should also exhibit the 

same behaviour. A first-order phase transition however is expected to become weaker 
under an increasing external field and it becomes a difficult task to recognize a 
double-peak in a weak first-order phase transition. Further increase in external field 
reduce the strength of the transition and the phase transition disappears for large 
external field strengths. This process is the result of two competing interactions: 
spins form sizable clusters in some regions because of the tendency to align in the field 
direction, while the thermal fiuctuations result in disintegration of the clusters. The 
external field yields the ordering of the system at higher temperatures, shifting the 
critical temperature to higher values as well as reducing fiuctuations in the system. 
On the other hand, the dynamical changes associated with a phase transition in the 
system are refiected in the fiuctuations and variations in the cluster size. The aim 
in this work is to observe, by using a cluster algorithm, the dynamical behaviour of 
cluster size variations with respect to temperature and the external field. 

In the present work, the critical behaviour of the two-dimensional Potts model with 
seven states under external magnetic field is studied using cluster algorithm. Equi- 
librium averages, fiuctuations in the average cluster size and the histograms for 
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energy, order parameter and the average cluster size are obtained as a function of 
the temperature for different field strengths. 



The Hamiltonian of the two-dimensional Potts model is given by 



(1) 



<i,j> i 



Here K = J / kT ; where k and T are the Boltzmann constant and the temperature 
respectively, and J is the magnetic interaction between spins cxj and aj, which can 
take values 0, 1, 2, g — 1 for the g-state Potts model and H = h/kT with h is the 
external field along the orientation 0. Reader can refer to the review article by Wu 
P for detailed information about the model. Order of the transitions can be studied 
by calculating specific heat 



and the Binder cumulant p 



on finite lattices, where E is the energy of the system. 

The dynamical evolution of clusters and fluctuations in the observables are strongly 
dependent on the correlations in the system. When the temperature is high the 
clusters are small, and they start growing as the critical point is approached. If 
the correlation length in the system is finite (especially when it is shorter than the 
lattice size), existing large clusters may break down to smaller ones with thermal 
fluctuations in the system. Hence in such a system, very large and very small clusters 
may coexist, resulting in large fluctuations in the cluster size. If the correlation 
length is very large as in the case of weak flrst-order phase transition or inflnite 
as in the case of a second-order transition, the existing large clusters can not easily 
disintegrate with thermal fluctuations. A similar effect can be seen when an external 
fleld is turned on in a system with flrst-order phase transition. The effect of the 
external magnetic fleld aligns the spins parallel to the fleld direction creating large 
clusters and preserve these clusters from the effects of thermal fluctuations. These 
considerations led us to the usage of the average cluster size, fluctuations in the 
cluster size, cluster size distribution as the operators to investigate the changes in 
the phase transition with respect to variations of temperature and external fleld 
strength. The cluster algorithm which introduces a global update by means of 
selecting clusters and updating all spins in the cluster proved to be very successful 
in eliminating critical slowing down and super critical slowing down as in the flrst- 
order transitions in systems such as Potts and 0(N) models. Studying any 
operator which uses information related to the clusters is extremely straightforward 
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since this information is inherited in the cluster algorithm. The algorithm used in 
this work is similar to Wolf's algorithm, with the exception that, H is incorporated 
following Dotsenko et al |rT[ and before calculating the observables, searching the 
clusters is continued until the total number of sites in all searched clusters is equal 
to or exceeds the total number of sites in the lattice. 

In this work, the two-dimensional, q = 7 Potts model with external magnetic 
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Fig. 1 . Energy histograms at the critical 
point Kc for several values of H. 



Fig. 2 . Specific heat peaks vs. H. 
Dashline is drawn to guide the eye. 



field has been simulated on 64 x 64 square lattice. After thermalization with 
10^ — 2 X 10^ sweeps, 5 x 10^ — 2 x 10^ iterations are performed at diff'erent 
values of the coupling K and the field H. Longer runs with up to 3 x 10^ it- 
erations are done near the finite-size critical value at each field strength H. 
For H = 0.00,0.01,0.015,0.016,0.02,0.03 and 0.05 the long runs are performed 
at = 1.2909, 1.2774, 1.2707, 1.2696, 1.2652, 1.2531 and 1.2299 respectively, where 
these points are chosen as the estimated peak positions of the specific heat. From 
these long runs we have evaluated by appropriate reweighting [12| the energy as a 
function of the temperature for each value of H. The continuous curves of energy, 
specific heat and the Binder cumulant obtained through extrapolation from one long 
run are in very good agreement with the data points over rather wide ranges of K 
for each H value . The point should be stressed here is that the success of an extrap- 
olation depends on number of statistically independent configurations used, which is 
provided by the cluster algorithm employed here. As a first criterion to distinguish 
phase transition, we have checked for each value of H, the energy distribution at 
the critical temperature where the specific heat possess a peak. At H = 0, where 
seven-state Potts model exhibits first-order character, a distinct double-peak is ob- 
served. As one can see from figure 1, with increasing H the double peak converges 
to a single peak while H reaches to the value of 0.016. Our simulations indicate that 
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H = 0.016 is almost at the verge of the first-order phase transition region and for 
larger values of H we have seen no double peak behaviour but from the observation 
of a single gaussian energy distribution, it is hard to conclude whether the transition 
becomes second-order or it is totaly wiped out. 
Prom specific heat peaks and the Binder cumulant minima, relevant information 
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Fig. 3 . Cluster size distributions for seve- 
ral values of H. 



Fig. 4 . Fluctuations in cluster size 
vs. K. 



for the phase transition may also be obtained. For a first-order phase transition, 
the specific heat exhibits Dirac-delta function hke shape. While the system moves 
towards the region of softer phase transitions with increasing iJ, the peak widens 
as well as the height is reduced. In the case of no phase transition, the system may 
still be expected to possess some fluctuations in the energy. For the set of values 
of the increasing external field, one obtaines a series of specific heat curves which 
become lower and wider as one moves towards the lower valus of K. Instead of a 
polynomial fit, here we have simply joined the tips of the successive specific heat 
curves by straight lines and displayed it in figure 2. This most naive presentation 
already prevails the existence of two distinct regimes and the turnover point is at 
B. — 0.016 where the line of first order transitions comes to an end. Our data on 
the Binder cumulant minima led to the same conclusion. 



In order to have more microscopic insight and demonstrate what possibly the cluster 
search may add to the knowledge, we studied the time evolution of average cluster 
size (CS), its fiuctuations (FSC) and the cluster size distributions (CSD). While the 
average energy fiuctuates within a limited range about a mean value, the size of 
then existing clusters may easily fluctuate between the two extremes, namely from 
a single spin to a cluster of lattice size. Hence it seems to us that one may draw 
more information from cluster formation about the dynamical changes occuring in 
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the system. Although the chistcrs wc consider are of the Swendsen-Wang type, since 
the changes in cluster size are due to the update algorithm used, the observables 
not dealing with the exact cluster size (hke FCS and CSD) should give correct 
information on the system. We have studied the time evolutions over 2-3 million 
iterations of average cluster sizes about Kf. for each value of the external field. 
These data are very much the same of the corresponding time-series for energy but 
more amplified and displays more insight about the grouping of spins leading to 
the concerned energy value. Two-state structures of first-order transitions and the 
long range fiuctuations of second order phase transitions are more easily detectable 
comparing to the energy-time series. In figure 3, we show the corresponding cluster 
size distributions in terms of histograms over the cluster size (normalized to the 
volume) in the range of to 1. From figure 3, one can easily follow the appearence 
of the two-state of small and large size clusters as for the H — 0.00 case, almost 
equal weight for all sizes as for H = 0.016 and loosing the small and large clusters in 
favor of medium sizes for H > 0.016. The point to bring into the readers attention 
here is that the distribution of small clusters (bins < 0.1) are very different for the 
three cases considered. For a first-order phase transition a substantial amount of 
small clusters are present and they gradually disappear while the phase transition 
weakens. 

The fiuctuations FCS {FCS =< {CSy > — < CS >^) in the average cluster size 
are calculated for different values of if , as a function of K. For the case of no 
external field, at small K values (high temperatures) the cluster sizes are small, 
hence FCS is small. As K increases, the cluster sizes are getting larger, but FCS 
increases due to the existence of the small-size clusters as well as the large ones. 
At the critical point, formation of the largest clusters leads to the largest value of 
FCS. One can make a similar discussion for K > Kc (when the critical point is 
approached from above). When the external field is turned on, spin alignments in 
the field direction reduces the fluctuations and the probability of finding large and 
small clusters in coexistence decreases with the increasing magnetic field. FCS for 
H = 0.00,0.01,0.016,0.02 and 0.05 are plotted in figure 4. Joining the tips of the 
successive FCS curves in figure 4 yields nothing but exactly what is displayed in 
figure 2 apart from a small shift in K due to the finite size effects. As can be seen from 
these plots, increasing magnetic field reduces the critical temperature, as well as the 
peak heights. As H increases, the FCS curves widen and after some H value, it is 
hard to see a distinct maximum point. At H = 0.05 the FCS curve is step function 
like where the phase transition is already wiped out. Besides the appearence of two 
different regimes on either side of the field value H — 0.016, figure 4 also displays 
differences in the distinct shapes of the FCS curves belonging to two regimes. We 
observed that the FCS curves for H > 0.016 are continuous and a polynomial fit 
is possible which may be denoted as the sign of a smooth change in the average 
cluster size in the system. The curves for H < 0.016 have the shape of a spike, for 
which the only possible polynomial fit is to have two distinct polynomials joining at 
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the tip. This appcarcnt discontinuity in the rate of fluctuations in cluster size may 
be attributed to the two state structure of first-order phase transitions. The low 
temperature asymptotic values for all H seem to be almost the same regardless of 
the order of transition. Because of small thermal fluctuations at low temperatures, 
the clusters freeze after they are formed, resulting in a small FCS value. 

In conclusion, one can see from CSD and FCS plots and from the energy histograms 
obtained near Kc for each considered value of H and from the extrapolations per- 
formed using the flttcd values of the maximums of the speciflc heat, the first-order 
transition in 2D g = 7 Potts model seems to disappear at H = 0.016. What we 
would like to bring into attention here is that the observation of the temperature 
dependence of fiuctuations in the average cluster size and the cluster size distribu- 
tions gives valuable information about the nature of the transition and this kind 
of investigation may enrich the physical insight when employed in studying phase 
transitions. 

Further work is planned to study random-field and random-bond Potts models using 
the same algorithm. 
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FIGURE CAPTIONS 

Figure 1. Energy histograms at the critical point Kc for several values of H. 
Figure 2. Specific heat peaks vs. H . Dashline is drawn to guide the eye. 
Figure 3. Cluster size distributions for several values of H. 
Figure 4. Fluctuations in cluster size vs. K. 
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